A numerical approach to spherical indentation techniques for material property evaluation

In this work, some inaccuracies and limitations of prior indentation theories, which are based on experimental observations and the deformation theory of plasticity, are investigated. Effects of major material properties on the indentation load-deflection curve are examined via finite element (FE) analyses based on incremental plasticity theory. It is confirmed that subindenter deformation and stress-strain distribution from deformation plasticity theory are quite dissimilar to those obtained from incremental plasticity theory. We suggest an optimal data acquisition location, where the strain gradient is the least and the effect of friction is negligible. A new numerical approach to indentation techniques is then proposed by examining the FE solutions at the optimal point. Numerical regressions of obtained data exhibit that the strain-hardening exponent and yield strain are the two key parameters which govern the subindenter deformation characteristics. The new indentation theory successfully provides a stress-strain curve and material properties with an average error of less than 3%.     

疲劳损伤伤演方程中材料常数的确定

损伤力学作为近年来固体力学的一个非常活跃的分支,其本构关系和损伤演化方程是核心内容。损伤演化方程中材料常数的确定,是损伤力学中最基本却很重要的工作。利用大范围损伤下的预估拉压和弯曲疲劳裂纹形成寿命的封闭解答,应用最小二乘法给出了确定损伤演化方程dD/dN=α(Δε)m中材料常数α、m的方法。该方法可作为确定这种类型损伤演化方程材料常数的一般方法。     

Elastic and Damage Longitudinal Shear Behavior of Highly Concentrated Long Fiber Composites

The elastic and damage longitudinal shear behavior of highly concentrated long fiber composites is analyzed by means of a simplified model where it is supposed that the fibers are rigid and touch each other in a regular hexagonal array. In the microscopic unit cell the problem is reduced to six similar problems of antiplane deformation on an equilateral circular triangle (see forthcoming Figure 2). These problems are solved in closed form by the complex variable method, and the solution is used to determine the longitudinal shear moduli, and to study their dependence on the microscopic damage caused by the circumferential debonding at the fiber–matrix interface. Subsequently, the damage evolution is investigated under the hypothesis that the microcracks propagate according to the Griffiths energy criterion. The elastic domain, where there is no damage propagation, is determined and it is shown that it is a polygonal convex set symmetric with respect to the origin. The overall damage evolution is discussed in detail and illustrated with some examples which highlight the very rich nature of the proposed model.     

Thermal buckling and post-buckling of pinned–fixed Euler–Bernoulli beams on an elastic foundation

In this article, both thermal buckling and post-buckling of pinned–fixed beams resting on an elastic foundation are investigated. Based on the accurate geometrically non-linear theory for Euler–Bernoulli beams, considering both linear and non-linear elastic foundation effects, governing equations for large static deformations of the beam subjected to uniform temperature rise are derived. Due to the large deformation of the beam, the constraint forces of elastic foundation in both longitudinal and transverse directions are taken into account. The boundary value problem for the non-linear ordinary differential equations is solved effectively by using the shooting method. Characteristic curves of critical buckling temperature versus elastic foundation stiffness parameter corresponding to the first, the second, and the third buckling mode shapes are plotted. From the numerical results it can be found that the buckling load-elastic foundation stiffness curves have no intersection when the value of linear foundation stiffness parameter is less than 3000, which is different from the behaviors of symmetrically supported (pinned–pinned and fixed–fixed) beams. As we expect that the non-linear foundation stiffness parameter has no sharp influence on the critical buckling temperature and it has a slight effect on the post-buckling temperature compared with the linear one.     

A note on the extensional viscosity of elastic liquids under strong flows

In this article a parametric study based on a balance between viscous drag and restoring Brownian forces is used in order to construct a nonlinear dumbbell model with a finite spring and a drag correction for a dilute polymer solution. The constitutive equations used are reasonable approximation for describing flows of very dilute polymer solutions such as those used in turbulent drag reduction. We investigate the response of an elastic liquid under extensional flows in order to explore the roles of a stress anisotropy and of elasticity in strong flows. It is found that for low Reynolds numbers, the extensional viscosity of a dilute polymer solution is governed by two parameters: a Deborah number representing the importance of the elasticity on the flow and the macromolecule extensibility that accounts for the viscous anisotropic effects caused by the macromolecule orientation. Two different asymptotic regimes are described.The first corresponds to an elastic limit in which the extensional viscosity is a function of the Deborah number and the particle volume fraction. The second is an anisotropic regime with the extensional viscosity independent of Deborah number but strongly dependent on macromolecule aspect ratio. The analysis may explain from a phenomenological point of view why few ppms of macromolecules of high molecule weight or a small volume fraction of long fibres produce important attenuation of the pressure drop in turbulent flows. On the basis of our analysis it is seen that the anisotropic limit of the extensional viscosity caused by extended polymers under strong flows should play a key role in the attenuation of flow instability and in the mechanism of drag reduction by polymer additives.     

Bounds on the Elastic Moduli of Completely Random Two-Dimensional Polycrystals

Explicit bounds on the elastic moduli of completely random planar polycrystals, the shape and crystalline orientations of the constituent grains of which are uncorrelated, are derived and calculated for a number of crystals of general two-dimensional anisotropy. The bounds on the elastic two-dimensional bulk modulus happen to coincide with the simple third order (in anisotropy contrast) bounds for the subclass of idealistic circular cell polycrystals. The bounds on the shear modulus are close to the much simpler bounds for circular cell polycrystals, which approximate aggregates of equiaxed grains.     

常速度拉伸实验惯性效应的实验研究

本文统计了１５０根Ａ３低碳钢比例拉伸试件的断口位置，并测定了常速度拉伸（４００ｍｍ／ｍｉｎ，０．５ｍｍ／ｍｉｎ）时应变随时间和沿杆长的分布以及应力应变关系；发现７２％试件的断口靠近不动端内。同一时刻在设定速度为４００ｍｍ／ｍｉｎ时杆中最大工程应变差１８％．     

Propagation of longitudinal elastic waves in composites with a random set of spherical inclusions (effective field approach)

The work is devoted to the problem of plane monochromatic longitudinal wave propagation through a homogeneous elastic medium with a random set of spherical inclusions. The effective field method and quasicrystalline approximation are used for the calculation of the phase velocity and attenuation factor of the mean (coherent) wave field in the composite. The hypotheses of the method reduce the diffraction problem for many inclusions to a diffraction problem for one inclusion and, finally, allow for the derivation of the dispersion equation for the wave vector of the mean wave field in the composite. This dispersion equation serves for all frequencies of the incident field, properties and volume concentrations of inclusions. The long and short wave asymptotics of the solution of the dispersion equation are found in closed analytical forms. Numerical solutions of this equation are constructed in a wide region of frequencies of the incident field that covers long, middle, and short wave regions of propagating waves. The phase velocities and attenuation factors of the mean wave field are calculated for various elastic properties, density, and volume concentrations of the inclusions. Comparisons of the predictions of the method with some experimental data are presented; possible errors of the method are indicated and discussed.     

On the application of the Wiener–Hopf technique to problems in dynamic elasticity

Many problems in linear elastodynamics, or dynamic fracture mechanics, can be reduced to Wiener–Hopf functional equations defined in a strip in a complex transform plane. Apart from a few special cases, the inherent coupling between shear and compressional body motions gives rise to coupled systems of equations, and so the resulting Wiener–Hopf kernels are of matrix form. The key step in the solution of a Wiener–Hopf equation, which is to decompose the kernel into a product of two factors with particular analyticity properties, can be accomplished explicitly for scalar kernels. However, apart from special matrices which yield commutative factorizations, no procedure has yet been devised to factorize exactly general matrix kernels.

This paper shall demonstrate, by way of example, that the Wiener–Hopf approximant matrix (WHAM) procedure for obtaining approximate factors of matrix kernels (recently introduced by the author in [SIAM J. Appl. Math. 57 (2) (1997) 541]) is applicable to the class of matrix kernels found in elasticity, and in particular to problems in QNDE. First, as a motivating example, the kernel arising in the model of diffraction of skew incident elastic waves on a semi-infinite crack in an isotropic elastic space is studied. This was first examined in a seminal work by Achenbach and Gautesen [J. Acoust. Soc. Am. 61 (2) (1977) 413] and here three methods are offered for deriving distinct non-commutative factorizations of the kernel. Second, the WHAM method is employed to factorize the matrix kernel arising in the problem of radiation into an elastic half-space with mixed boundary conditions on its face. Third, brief mention is made of kernel factorization related to the problems of flexural wave diffraction by a crack in a thin (Mindlin) plate, and body wave scattering by an interfacial crack.     

Constitutive framework optimized for myocardium and other high-strain, laminar materials with one fiber family

Central to this analysis is the identification of six rotation invariant scalars α_1-6 that succinctly define the strain in materials that have one family of parallel fibers arranged in laminae. These scalars were chosen so as to minimize covariance amongst the response terms in the hyperelastic limit, and they are termed strain attributes because it is necessary to distinguish them from strain invariants. The Cauchy stress t is expressed as the sum of six response terms, almost all of which are mutually orthogonal for finite strain (i.e. 14 of the 15 inner products vanish). For small deformations, the response terms are entirely orthogonal (i.e. all 15 inner products vanish). A response term is the product of a response function with its associated kinematic tensor. Each response function is a scalar partial derivative of the strain energy W with respect to a strain attribute. Applications for this theory presently include myocardium (heart muscle) which is often modeled as having muscle fibers arranged in sheets. Utility for experimental identification of strain energy functions is demonstrated by showing that common tests on incompressible materials can directly determine terms in W. Since the described set of strain attributes reduces the covariance amongst response terms, this approach may enhance the speed and precision of inverse finite element methods.